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Abstract 

The hyperfine splitting of the ground state of H-, Li-, and B-like ions is investigated in details within 
the range of nuclear numbers Z = 7 — 28. The rigorous QED approach together with the large-scale 
configuration-interaction Dirac-Fock-Sturm method are employed for the evaluation of the interelectronic- 
interaction contributions of first and higher orders in l/Z. The screened QED corrections are evaluated to all 
orders in aZ utilizing an effective potential approach. The influence of nuclear magnetization distribution is 
taken into account within the single-particle nuclear model. The specific differences between the hyperfine- 
structure level shifts of H- and Li-like ions, where the uncertainties associated with the nuclear structure 
corrections are significantly reduced, are also calculated. 

PACS numbers: 32.fO.Fn, 31.15.aj, 31.30.J- 
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I. INTRODUCTION 



Accurate knowledge of the hyperfine structure lines of middle- Z multicharged ions is of great 
interest due to suggested observations of these lines from hot rarefied astrophysical plasmas yj,|2fl. 
Such observations may allow one to study the chemical and isotopic compositions of the supernova 
remnants, and the hot interstellar medium, including the galactic halos, which are the main types 
of objects from which intense emission lines are expected. The experiments on the determination 
of hyperfine splittings will also enable us to refine the deduction of nuclear magnetic moments of 
different isotopes and to inspect the various computational models employed for the theoretical de- 
scription of nuclear effects. High-precision measurements of theground-state hyperfine structure 
of heavy highly charged ions have been performed in Refs. [3, . Extension of these exper- 

iments to Li-like ions presently being prepared will provide tests of quantum electrodynamics 
(QED) in strong electric and magnetic fields on level of a few percent in specific difference of the 
hyperfine splitting values of H- and Li-like ions Q9Q. In this difference the main theoretical uncer- 
tainty which originates from the nuclear magnetization distribution correction (Bohr-Weisskopf 
effect) is essentially reduced. In specific differences of heavy H- and B-like ions or Li- and B-like 
ions the same reduction of the theoretical uncertainty can be also achieved. This becomes clear 
from the approximate analytical expressions for the Bohr-Weisskopf correction given in Ref. 11 1 CXI . 

The theoretical investigations of the hyperfine splitting of H- and Li-like multicharged ions 
in the middle-Z region have some history. The first accurate calculation (~ 0.1%), based on a 
combination of 1 jZ perturbation theory and the nonrelativistic configuration- interaction Hartree- 
Fock method, was performed in Refs. lllll.ll2n. Later, Boucard and Indelicato fl 1 3M employing the 
multi-configuration Dirac-Fock method presented the evaluation of the hyperfine splitting values 
over the entire range of the nuclear charge numbers Z = 3 — 92. Expansion in aZ of the QED 
correction has been worked out in Refs. Jl^lfisi lfll (for earlier studies see references therein and 



recent reviews 1117 . 



180). However, the application of the aZ-expansion is restricted to s-states 



in one-electron ions and limited by its convergence property. Therefore, we evaluate the radiative 
corrections numerically to all orders 'maZ accounting for the interelectronic-interaction effects by 
means of local screening potentials. All-order calculations of one-loop QED contributions to the 
hyperfine structure for middle- Z ions have been previously performed for the Is state [Jjl 



22fl, for the 2s state Q21 



20, 



21, 



2211 . and for the 2pi/ 2 state [23]. However, almost all these calculations 



of the QED corrections were dealing with one-electron ions only, where the screening effects are 
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absent. 



In the present paper, we calculate the ground-state hyperfine structure of H-, Li-, and B-like 
sequences in the middle-Z region. The one-loop radiative corrections are evaluated to all orders 
in olZ employing an effective local screening potential. Many-body effects are taken into account 
to the first order in 1/Z within the QED perturbation theory and to higher orders within the large- 
scale configuration-interaction Dirac-Fock-Sturm method (CI-DFS). The single-particle nuclear 
model is employed for the evaluation of the Bohr-Weisskopf correction. The main goal of this 
work is to improve the accuracy of previous results for the hyperfine structure of H- and Li-like 
ions and to present novel calculations for the B-like sequence in the middle-Z region. 

The paper is organized as follows: In the next section the basic formulas for the hyperfine split- 
ting are given and the derivation of the various contributions is described. In Section [Til] we present 
the numerical results for all contributions and compare the total values with previously reported 
calculations and with existing experimental data. Section [IV] provides a complete compilation of 
the total values for the hyperfine splitting of H-, Li, and B-like ions as well as the results for the 
specific differences between the hyperfine structure of H- and Li-like ions. We close with a short 
summary and point out the main achievements of the present work. 

Relativistic units (h — 1, c = 1, m — 1) and the Heaviside charge unit [a = e 2 / (47r), e < 0] 
are used throughout the paper. 

II. BASIC EXPRESSIONS 

The interaction of atomic electrons with the nuclear magnetic-dipole moment is described by 
the Fermi-Breit operator, which is conveniently written as a scalar product of two tensor operators 



where /x is the nuclear magnetic moment operator acting in the space of nuclear states. The 
electron part T is defined by the following expression 



where index i refers to the i-th electron of the atom, a is the Dirac-matrix vector, and rij = r^/r^. 
This interaction leads to the hyperfine splitting of the atomic levels. For an ion with one electron 




(1) 




(2) 
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(e.g., ns or npi/ 2 state) over the closed shells this splitting can be written in the form 

AE (a) = <*{<*Z)*9l 21+1 1 



A{otZ){l -8){l-e) + ^B(aZ) + -^C(Z, aZ) + x rad 



mp ( J + l)(2/ + l)(l + f)3 

1 1 

(3) 

Here Z is the nuclear charge number, m v and M are the proton and nuclear masses, respectively. 
Within the approximation of noninteracting electrons, where the contribution of the closed shells 
is neglected, the hyperfine splitting is explicitely determined by the quantum numbers of valence 
electron state a, which is characterized by the principal quantum number n, the angular momentum 
j, its projection rrij, and the parity I. A nucleus with spin / possesses a nuclear g factor gj = 
fi/ [xnL where fi is the nuclear magnetic moment and is the nuclear magneton. A(aZ) is the 
one-electron relativistic factor, 5 and e are, respectively, the corrections for distributions of the 
charge and magnetic moment over the nucleus; the functions B(aZ) and C(Z, aZ) determine the 
corrections for the electron-electron interaction of first and higher orders in l/Z, respectively; x ra( j 
is the QED correction. These terms are subsequently described in the following subsections. 

A. One-electron contributions 



The relativistic factor A(aZ) corresponding to the point-like nucleus is known analytically [12411 



n 3 (2l + l)42K( 1 + n r )-N} 
A{aZ) = iV- 7 (4 7 2 - 1) ' (4) 

where n r = n — \k\ is the radial quantum number, k = (—iy +l+1 ^ 2 (j + 1/2), 7 = a/k 2 — (aZ) 2 , 
N = ^n? + 2n r 7 + k 2 . The nuclear charge distribution correction 5 can be found either ana- 
lytically JloL I25I1 or numerically by solving the Dirac equation with the Coulomb potential of the 
extended nucleus. In this work it is evaluated numerically employing the homogeneously-charged- 
sphere model for the nuclear charge distribution. In order to estimate the uncertainty due to the 
model dependence the Fermi model is used as well. The Bohr-Weisskopf correction e originates 
from the spatial distribution of the magnetic moment inside the nucleus. For a rigorous treatment 



of this effect for low-Z systems we refer to Ref. [26]. In the present work we restrict our consider 



ation to models in which it can be accounted for by replacing the factor 1/r 2 in Eq. © by F(r)/r 2 , 
where F(r) is the volume distribution function. For example, in case of the sphere model it reads 

F(r)=< (£) ' r - fl ° , (5) 
1, r > R 



where R = a/5/3 (r 2 ) 1 / 2 is the radius of the sphere, and (r 2 ) 1 / 2 is the charge root-mean- square 
radius of the nucleus. However, with the sphere model one can not always describe adequately 
the nuclear magnetization distribution. The approximation of the nuclear single-particle model 
is widely used for the evaluation of the Bohr-Weisskopf correction 

Within this model the nuclear magnetization is determined by the total angular momentum of the 
unpaired nucleon (proton or neutron). Accordingly, the nuclear g factor gj is just the Lande factor 
of an extra nucleon, which is defined by the well-known formula 



gi = fi/ fi N I 



{9l + 9s) + (9l ~ 9s) 



L(L + l)-3/4 



(6) 



where L is the nuclear orbital momentum, g^ and g$ are the orbital and spin g factors of the valence 
nucleon, respectively. In case of a valence proton gi = 1, while for an extra neutron gi = 0; (75 is 
chosen such as to reproduce the experimental value of the nuclear magnetic moment fi according 
to Eq. ©. For nuclei with odd or even nuclear charge numbers the role of the unpaired nucleon 
is either played by a proton or a neutron, respectively. In the framework of the nuclear single 



29 



particle model the radially symmetric distribution function F(r) has been derived in Refs. [27, 
301. Here we neglect the contribution of the spin-orbit interaction and employ the homogeneous 

nn 

distribution for the radial part of the odd nucleon wavefunction inside the nucleus [10, 30]. In this 
approximation F(r) reads 



F(r) 
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8(1 + 1 



:9s 



1 



9l 



r<R , (7) 



27 + 3 7(27 + 3) 
;9s + " T , ' ql 



8(7 + i; 



2(7+1) 



r < Ro 



(8) 



for 7 = L — |. For r > R the distribution function F(r) = 1. In the case of 14 N with 7 = 1 we 
follow the work 111 ill and assume that the nuclear magnetization is determined by the odd proton 
and neutron. The corresponding formulas for e were derived in Refs. ft 1 lL 13111 . The uncertainty 
of the Bohr-Weisskopf correction is estimated as the maximum of two values: 50% of e itself 
and the difference between e obtained in single-particle and sphere nuclear models. As in our 
previous studies (see, e.g., the related discussion in Ref. |27|]), the uncertainty obtained by this 
procedure must generally be considered only as the order of magnitude of the expected error bar. 
More accurate calculations of the Bohr-Weisskopf effect must be based on many-particle nuclear 
models and should include a more rigorous procedure for determination of the uncertainty. The 



nuclear vector polarizability correction derived in Ref. Il26ll is assumed to contribute less than the 
uncertainty of e indicated above. 

Separating out the nuclear parameters and the nonrelativistic value of the hyperfine splitting 
one finds that the one-electron contributions considered above can be numerically evaluated in 
terms of a matrix element 

A(aZ)(l-5)(l-e) = G a (a\T \a) (9) 

of the zero component T of the operator T given by Eq. 0, multiplied by the magnetization 
distribution function F(r). The wavefunction \a) of the valence state, characterized by quantum 
numbers a = n, j, rrij, and I, is obtained as a solution of the Dirac equation with the potential of 
the extended nucleus. The multiplicative factor G a reads 

_ n 3 (2Z + l)j(j + l) 
Ga ~ 2(aZfm 3 • (10) 

B. Many-electron contributions 

Now we pass to the many-electron corrections. The term B(aZ)/Z in Eq. ([3]) determines 
the interelectronic-interaction correction of the first order in 1/Z. A rigorous QED treatment 



of this contribution can be carried out utilizing the two-time Green's function method [32]. To 
simplify the derivation of formal expressions, it is convenient to incorporate the core electrons as 
belonging to a redefined vacuum. This leads to merging the interelectronic-interaction correction 



of order 



Refs. [33 



/Z with the one-loop radiative corrections. Such a treatment was applied previously in 



3411 . The corresponding expression for the interelectronic-interaction correction reads 



B(aZ)/Z = 2G a J2 E 



a (ac\I(0)\nc)(n\To\a) e ^t^ c (ac\I(0)\an}(n\T \c) 



_ £ \p a (ac\I(e a - e c )\cn)(n\T \a} _ (ac\I(e a - e c )\na)(n\T \c) 

n n 

~ [(a\T \a) - (c\T \c)) (ac\I'(e a - e e )\ca) \ , (11) 

where e m are the one-electron energies, J (a;) = e 2 a^a u D^uj), = dl(u)/duj, = (1, a), 
and D^ u (uj) is the photon propagator. It should be noted that the total 1/Z interelectronic- 
interaction correction given by Eq. ([111) is gauge independent. We perform the calculation employ- 



6 



ing Coulomb and Feynman gauges for the photon propagator, thus receiving an accurate check of 
the gauge invariance of the results. 

The interelectronic-interaction correction of higher orders C(Z,aZ)/Z 2 is calculated within 
the framework of the large-scale configuration-interaction method in the basis of Dirac-Fock- 
Sturm orbitals fl35ll . This method was successfully employed in our previous atomic calculations 

2 



[3 



37, 



38, 



39 



The interelectronic-interaction operator employed in the Dirac-Coulomb-Breit 



equation reads 



Vw 



Cti-Ctj (dti • Ti)(0Lj ■ Tj] 



2r, ; 



2r 3 - 



(12) 



where the sum runs over all electrons. A scaling parameter A is introduced to separate terms of 
different order in 1/Z from the numerical results with different A. This representation allows us to 
perform the expansion in powers of A. In this way, the higher-order term is written as 



C{Z,aZ)/Z 2 = G a { (^ A ( 7 JM J )|To|^ A ( 7 JM J )) - (* A ( 7 JMj) |r |* A ( 7 JMj)> 



- — (^xilJMj^T^^JMj)) 



A=l 



A=0 



A=0 



(13) 



The many-electron wavefunction \T / A ( 7 JMj) is characterized by the total angular momentum J, its 
projection Mj, and the rest quantum numbers j. The configuration-interaction matrix contains all 
single, double, and triple positive-energy excitations. Single-electron excitations to the negative- 
energy spectrum were accounted for in the many-electron wavefunction \I/ A ( 7 JMj) employing 
perturbation theory. 

The calculation of the interelectronic-interaction corrections B(aZ)/Z and C(Z, aZ)/Z 2 is 
performed employing the homogeneously-charged-sphere model for the nuclear charge distribu- 
tion and single-particle model for the nuclear magnetic moment distribution. 



C. One-loop radiative contribution 

The one-loop radiative contribution x ra d appears as the sum of vacuum-polarization (VP) and 
self-energy (SE) corrections, x rad = x V p + x SE , as depicted diagrammatically in Figs. Q] and [2l re- 
spectively. However, for the Li- and B-like ions along with the one-electron part the correction x ra( j 
contains also the many-electron part. In order to account for many-electron effects we consider 
an effective spherically symmetric potential V cS that partly takes into account the interelectronic 
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FIG. 1: Feynman diagrams representing the vacuum-polarization correction to the hyperrine splitting. The 
wavy line indicates the photon propagator and the double line indicates the bound-electron wavefunctions 
and propagators. The dashed line terminated with the triangle denotes the hyperrine interaction. 
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FIG. 2: Feynman diagrams representing the self-energy correction to the hyperfine splitting. Notations are 
the same as in Fig.Q] 

interaction between the valence electron a and the core electrons c. This can be achieved by means 
of the Kohn-Sham screening potential derived within the density-functional theory 114 IN 



1 2 a 

V eS (r) = Kuc(r) + a I dr'—p t (r') - 



3r V32tt 2 



rp t [r / 



1/3 



(14) 



m, 



44 



450. Here V nuc is the po- 



which we employed successfully in previous calculations 11421 
tential of the extended nucleus and p t denotes the total one-electron density. In order to estimate 
the sensitivity of the result on the specific choice of the screening potential we consider also the 
core-Hartree potential. 

The VP correction x V p is divided into the electric-loop part, Fig. [T](a,b), which accounts for 
the VP correction to the scalar binding potential V cS , and the magnetic-loop part, Fig.[T](c), corre- 
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sponding to the VP-corrected hyperfine interaction potential. The expression for the electric-loop 
term Fig.[T](a,b) reads 



2G a 



(a\T \n) (n\Uy P \a) 



(15) 



where Uyp represents the renormalized one-loop VP potential. It is divided into the Uehling and 
Wichmann-Kroll parts, U$ P = U^' el + L^" 61 . The Uehling part can be evaluated according to 
the well-known equation 

2a 2 Z 



U. 



Ue-el, 
VP I 



37T 



-2 r-r' \t 



(16) 



where the density p e ff is related to the effective binding potential V c f[ [ via the Poisson equation 
AVrfr(r) = AnaZ p e fr(r) ]. The Wichmann-Kroll part can be generated by summing up the partial- 
wave differences between the unrenormalized total VP potential and the unrenormalized Uehling 



term 



470. hi this work we emp loy the approximate formula for the Wichmann-Kroll electric - 



loop potential derived in Ref. |48|]. The correction to the hyperfine splitting due to the magnetic 
loop x™p can be written in the form 



Xyp — G a (a\Uy-p\a) , 



(17) 



where f/yp is the VP-corrected hyperfine potential T . It can be renormalized utilizing the same 
scheme as for the electric loop. For the distribution function F(r) corresponding to the sphere 
model (see Eq. ©) we obtain the following analytical expression for the magnetic-loop Uehling 
term 



Ue— ml / 
VP \ 



a [n x a] 3 



7T 



16iS 



ArRc 



P 1 (R + r)+ i3 1 {\R - r\) + 2(R + r) f3 2 {R + r 



-2|i2o - r\ fc(\Ro - r\) + (3 3 (R + r) - p 3 (\R - r\) 



(18) 



where the function (3 n is defined as 

c\ poo 

Pn(r) = ~ 



dt 



t n+2 y L + 2ty 



(19) 



The contribution of the remaining Wichmann-Kroll magnetic-loop term is relatively small. For 
the case of the Is and 2s states the values for the Wichmann-Kroll magnetic-loop term are taken 
from Ref. Q22Q . For 2s we have also incorporated the screening effect, assuming that the screening 



coefficient is the same as for the Uehling magnetic-loop term. For the 2p 1 / 2 state Xy^~ ml turns 
out to be smaller than the uncertainties assigned to the calculation. 

Now let us turn to the evaluation of the SE correction. The formal expression can be derived by 
means of the two-time Green's function method Q32Q. The SE correction appears as the sum, xse = 
x se + x se + X SE' °f me irreducible Xg E , reducible , and vertex Xg|^ terms, respectively. The 
irreducible part, depicted in Fig.[2](a,b) with the intermediate state energy e n ^ e a , is represented 
by the expression 



x 



SE 



en ^ Sa (a\To\n)(n\ S(e a ) — 7°5m |a) 
2G a " — 



(20) 



where 5m is the mass counter-term and S(e) denotes the unrenormalized self-energy operator with 
matrix elements defined by 



{a\E(e)\b) 



% 

2^ 



(an\I (uj)\nb) 



e — uj — e n (l — iO) 



(21) 



Accordingly, the irreducible contribution can be written as nondiagonal matrix element of the 



self-energy operator. Thus, the renormalization scheme deve 



correction can be also applied in this case (see, e.g., Refs. |49|, 



oped for the first-order self-energy 



50, 



5111 ). Only a slight extension of 



the corresponding formulas for the case of a nondiagonal matrix element is needed. 
The expression for the reducible term is given by 

,dE(e) 



x 



„red 
SE 



G a (a\ 



de 



\a) (a|To|a) 



(22) 



while the vertex part, Fig.|2](c), reads 

x ™ = G a — / duV — 



(an 2 \I(uj)\n 1 a} (ni\T Q \n 2 ) 



m n 2 



u - e ni (l - i0))(e a - uj - e na (l - iO)) 



(23) 



Both reducible and vertex terms are ultraviolet-divergent. In order to isolate the divergencies in 
a covariant way, we separate out the zero-potential terms in which bound electron propagators 
are replaced by free propagators. The sum of the latter terms for the reducible and the vertex 
part is denoted by x^^ = x T ^ ^ + £gg • Their evaluation is performed in momentum space, 
where the ultraviolet divergencies can be canceled in a standard way. The remaining part of the 
reducible and the vertex contribution Xg^ 1 ^ is ultraviolet finite. However, we note that the term 



with e 



Hi 



- )(2 



e a in Eq. (1231) involves an infrared divergency, which is canceled by the cor- 



responding term of the reducible contribution. Performing the integration over the energy of the 



10 



virtual photon in these terms analytically we explicitly achieve the finite result in the sum of the 
reducible and vertex contributions. The remaining part of the many-potential term Xg^ 1+ ^ is in- 
frared finite, and can be calculated in coordinate space by means of a point-by-point subtraction 
of the corresponding contributions with free propagators inside of the self-energy loop. Angular 
integration and summation over intermediate angular momentum projections is carried out in a 



standard way H5 111 . The evaluation of the many-potential terms performed in the coordinate space 
involves an infinite summation over the angular-momentum quantum number k of intermediate 
states. This sum was extended up to |/t max | = 10 and the remaining part of the sum is estimated 
by a least-square inverse-polynomial fitting. One also observes that the results of the radial inte- 
gration converge better, when an extended model for the magnetization distribution is employed. 
In our calculations of the radiative corrections we have utilized the sphere model for the magnetic 
moment distribution function F(r). 

IE. NUMERICAL RESULTS 

Now let us pass to the presentation of the numerical procedure and the results for H-, Li-, 
and B-like sequences. The infinite summations over the complete spectrum of the Dirac equation 
involved in the numerical evaluations are performed employing the finite-basis set approach. The 
B-splines basis set was constructed utilizing the dual kinetic balance approach p2Q. The latter 
treats large and small components on equal footing and respects the charge conjugation symmetry. 
As a consequence no unphysical spurious states appear and moreover, it improves the convergence 
properties and the accuracy considerably. The values of the nuclear root-mean-square radii are 



taken from the tabulation 115 311 . The root-mean- square radius, in particular, for the 33 S isotope is 
assumed to be the average of the values given for even isotopes 32 S and 34 S, respectively. Empirical 
data for the nuclear properties: spin /, parity n, and magnetic moment /x//ijv are taken from 
Ref. [54]. All these values are also compiled in Table U We indicate the uncertainties assigned 
to the nuclear magnetic moments only if they exceed the level of 10~ 5 in the relative units. One 
has to note here, that the magnetic moment values obtained via the nuclear magnetic resonance 



technique do not usually account for the chemical shift 115 511 . which is of the order 10 3 — 10 4 or 
sometimes even larger. 
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A. H-like ions 



The individual contributions to the hyperfine splitting for the light H-like ions are presented 
in Table HI The values obtained for the radiative corrections ^se and xyp are in good agreement 
with the most accurate nonperturbative results Il22ll based on the Coulomb-Dirac Green function. 
The slight difference is explained by the finite-nuclear- size effects accounted for in the present 
work. The latter is especially important for the evaluation of the specific difference between H- 
and Li-like hyperfine splitting values, where the QED corrections have to be calculated within 
the same nuclear model. As one can see from the table the main uncertainty originates from the 
Bohr-Weisskopf correction e. In TableHUthe predictions for the total transition energies AE^ are 
compared with the results of previous calculations lllll. 1 1 3M . Deviations between our results and 
those reported in Ref. [| 1 3fl arise from the different treatment of the Bohr-Weisskopf effect. In work 
fl 1 311 the simple spherical model for the magnetization distribution was employed. Theoretical 
values for the total transition energies AE^ and wavelengths \^ are presented in Table IXl 
below. 



B. Li-like ions 



The consideration of Li-like ions we start with the results for the screened radiative corrections 
xse and xyp. The one-loop QED correction is conveniently represented in terms of the function 
D Tad defined as 



a n 

7T 



(24) 



In Tables IXTT1 and HVl the numerical results for individual contributions to the self-energy (-Dse) and 
vacuum-polarization (Dy P ) corrections are presented, respectively, for Z = 10, 15, 20, 25. 

Table IVl displays the individual contributions to the hyperfine splitting of the light Li-like ions. 
As in the case of H-like ions the main uncertainty originates from the Bohr-Weisskopf correction e. 



13|] account for the ra- 



Earlier calculations on the hyperfine structure of light Li-like ions [|11L 
diative correction on the basis of analytical expansion with respect to aZ. Here we have performed 
exact (to all orders in aZ) evaluations of one-loop QED corrections with an effective screening 
potential and with a nuclear vector potential involving an extended magnetization distribution. As 



compared to the results of works lllll, Il2h . several additional improvements have been achieved: 
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TABLE I: Individual contributions to the ground-state hyperfine splitting of the hydrogenlike ions. 
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I 71 
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1+ 


0.40376 


2.5579 


1.00393 


0.00067(3) 


-0.00004(28) 


0.00009 


0.00028 




\ 

2~ 


-0.28319 


2.6061 


1.00393 


0.00068(3) 


0.00114(88) 


0.00009 


0.00028 


17q7+ 


t; 

-4 


-1.8938(1) 


2.6953 


1.00514 


0.00081(3) 


0.00033(17) 


-0.00006 


0.00032 


19pS+ 


i 

5+ 


2.6289 


2.8976 


1.00651 


0.00099(3) 


0.00036(18) 


-0.00022 


0.00036 


21 Ne 9+ 


2 T 


-0.66180(1) 


2.9672 


1.00805 


0.00113(4) 


0.00058(29) 


-0.00037 


0.00040 


23 Na iO+ 


1 + 


2.2175 


2.9936 


1.00975 


0.00127(4) 


0.00035(18) 


-0.00053 


0.00044 


25 Mg ll + 


1 + 


-0.85545(8) 


3.0280 


1.01163 


0.00141(4) 


0.00057(29) 


-0.00068 


0.00049 


27^jl2 + 


- + 

2 r 


3.6415 


3.0605 


1.01367 


0.00156(5) 


0.00048(24) 


-0.00084 


0.00053 


29 g|134 


1 

5+ 


-0.55529(3) 


3.1168 


1.01589 


0.00173(5) 


0.00063(31) 


-0.00099 


0.00058 


31 pl44 
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1 + 


1.1316 


3.1888 


1.01828 


0.00191(5) 


0.00069(35) 


-0.00115 


0.00062 
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— 4- 
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0.64382 


3.2727 


1.02085 


0.00212(6) 


0.00106(53) 


-0.00130 


0.00067 


OK 1 (31 

OJ pi IDt 


t , 
— 4 

2 r 


0.82187 


3.3652 


1.02360 


0.00234(6) 


-0.00026(110) 


-0.00146 


0.00071 


37q16+ 


Q 

-4 

2 ^ 


0.68412 


3.3840 


1.02360 


0.00235(6) 


-0.00055(140) 


-0.00146 


0.00071 


39 


Q 

-4 

2 ^ 


0.39147 


3.4346 


1.02964 


0.00274(6) 


-0.0021(31) 


-0.00177 


0.00081 


41 j^18+ 


1+ 


0.21487 


3.4514 


1.02964 


0.00275(6) 


-0.0050(60) 


-0.00177 


0.00081 


43 Ca 19 + 


2 ~ 


-1.3176 


3.4928 


1.03294 


0.00297(7) 


0.00119(60) 


-0.00193 


0.00086 


45 Sc 20 + 


7 

2 ~ 


4.7565 


3.5443 


1.03644 


0.00320(7) 


0.00092(46) 


-0.00208 


0.00091 


47 Yj21+ 


5 

2 ~ 


-0.78848(1) 


3.5944 


1.04012 


0.00345(7) 


0.00160(80) 


-0.00224 


0.00096 


49^21+ 


7 

2 _ 


-1.1042 


3.5735 


1.04012 


0.00343(7) 


0.00137(69) 


-0.00224 


0.00096 


51y22+ 


7 

2 


5.1487 


3.5994 


1.04401 


0.00367(8) 


0.00107(54) 


-0.00240 


0.00101 


53Q.23+ 


3 
2 


-0.47454(3) 


3.6588 


1.04810 


0.00395(8) 


0.00149(75) 


-0.00256 


0.00106 


55 Mn 24+ 


5 
2 


3.4687 


3.7057 


1.05239 


0.00423(8) 


0.00109(54) 


-0.00272 


0.00111 


57p e 25+ 


1 
2 


0.090623 


3.7534 


1.05689 


0.00453(9) 


0.00279(140) 


-0.00289 


0.00117 


59 Co 26+ 


7 
2 


4.627(9) 


3.7875 


1.06161 


0.00483(9) 


0.00133(66) 


-0.00305 


0.00123 


61 Ni 27+ 


3 
2 


-0.75002(4) 


3.8221 


1.06655 


0.00514(9) 


0.00191(95) 


-0.00322 


0.00128 



TABLE II: Comparison of the ground-state hyperfine splitting of H-like ions between different 

theoretical calculations. The values of the energies are given in meV. 



Ion 




this work 


ni] 


[13J 


14 N 6+ 


0.40376 


0.21936(6) 


0.21937(4) 


0.21931 


27^12+ 


3.6415 


10.215(2) 




10.215 


45 Sc 20+ 


4.7565 


54.613(25) 




54.601 


57p e 25+ 


0.090623 


3.5109(49) 




3.5152 



the interelectronic-interaction corrections B(aZ)/Z and C(Z, aZ)/Z 2 have been calculated tak- 
ing into account explicitly the extended nuclear charge and magnetic moment distribution effects, 
moreover, the term C (Z, aZ)/Z 2 has now been evaluated within the framework of the relativistic 
CI-DFS method. The latter is especially important for ions of the higher Z re gion . The remaining 



one-electron corrections A(aZ), S, and e coincide with the results of works ft 1 X . 



1211 . The recoil 



effect for Li-like ions is partly accounted for by a factor (1 + m/M) 3 in Eq. ([3]). However, addi- 
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TABLE III: Individual contributions to the screened self-energy correction for the ground-state hyperfine 
structure of the lithiumlike ions, in units of the function Z?se- 



10 


-0.161 


2.723 


-2.661 


-0.100 


15 


-0.296 


2.685 


-2.760 


-0.372 


20 


-0.443 


2.543 


-2.772 


-0.672 


25 


-0.602 


2.380 


-2.775 


-0.999 



TABLE IV: Individual contributions to the screened vacuum-polarization correction for the ground-state 
hyperfine structure of the lithiumlike ions, in units of the function -Dyp. 



Z D 



Ue-el 
VP 



n Ue-ml n WK-el 



D 



WK-ml 
VP 



D 



VP 



10 


0.067 


0.058 


-0.000037 


-0.00020 


0.125 


15 


0.119 


0.098 


-0.00015 


-0.00076 


0.216 


20 


0.178 


0.140 


-0.00039 


-0.0020 


0.316 


25 


0.247 


0.186 


-0.00083 


-0.0041 


0.428 



tional contributions arising from the specific mass shift and spin-orbit recoil corrections [56] are 
significantly smaller than the uncertainty assigned for the Bohr-Weisskopf correction. 

In Table [VI] results for the total ground-state hyperfine splitting values of lithiumlike ions 
AE^ of different theoretical calculations are compared. In addition the experimental value of 
a recent measurement of the hyperfine splitting of lithiumlike 45 Sc 18+ ion, performed by resolv- 
ing the 2s hyperfine structure in the dielectronic recombination spectrum Il57ll . is given as well. 
The deviations between our results and values reported in Ref. [|13|l are mainly determined by 
interelectronic-interaction effects. The total theoretical values of the energies AE^ and wave- 
lengths A^ 2s ) are reported in Table IXl 



C. B-like ions 



Let us now turn to boronlike ions. The corresponding screened radiative corrections xse and 
xyp expressed in terms of the functions D SE and Dy P defined by Eq. (|24|) are presented in Ta- 
bles IVIII and IVIIIl respectively. The uncalculated Wichmann-Kroll contribution of the magnetic- 
loop D^^~ ml is of the same order as the corresponding electric-loop term D^^~ el . As can be 
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TABLE V: Individual contributions to the ground-state hyperfine splitting of the lithiumlike ions. The values 
of the nuclear parameters are the same as in Table HI 



Ion 


A(aZ) 


5 


e 




Xyp 


B(aZ)/Z 


C(Z,aZ)/Z 2 


14 N 4+ 


1.00557 


0.00067(3) 


-0.00004(28) 


0.00008 


0.00017 


-0.38146 


0.01800 


15^4+ 


1.00557 


0.00068(3) 


0.00114(88) 


0.00008 


0.00017 


-0.38101 


0.01798 


17 5+ 


1.00728 


0.00081(3) 


0.00033(17) 


-0.00001 


0.00021 


-0.33424 


0.01387 


19p6+ 


1.00923 


0.00099(3) 


0.00036(18) 


-0.00012 


0.00025 


-0.29767 


0.01104 


21 Ne 7+ 


1.01142 


0.00113(4) 


0.00058(29) 


-0.00023 


0.00029 


-0.26845 


0.00900 


23 Na 8+ 


1.01384 


0.00127(4) 


0.00035(18) 


-0.00035 


0.00033 


-0.24471 


0.00749 


25 Mg 9+ 


1.01650 


0.00141(4) 


0.00057(29) 


-0.00047 


0.00037 


-0.22488 


0.00634 


27 A1 10 + 


1.01941 


0.00156(5) 


0.00048(24) 


-0.00060 


0.00042 


-0.20823 


0.00544 


29 si ll+ 


1.02257 


0.00173(5) 


0.00063(32) 


-0.00073 


0.00046 


-0.19395 


0.00473 


31pl2+ 


1.02597 


0.00192(5) 


0.00070(35) 


-0.00086 


0.00050 


-0.18164 


0.00415 


33<jl3+ 


1.02963 


0.00212(6) 


0.00107(53) 


-0.00100 


0.00055 


-0.17086 


0.00368 


35q14+ 


1.03355 


0.00235(6) 


-0.00026(111) 


-0.00113 


0.00059 


-0.16166 


0.00329 


37 C1 14+ 


1.03355 


0.00236(6) 


-0.00055(140) 


-0.00113 


0.00059 


-0.16171 


0.00329 


39 K 16+ 


1.04219 


0.00275(6) 


-0.0021(31) 


-0.00142 


0.00069 


-0.14619 


0.00268 


41 K 16+ 


1.04219 


0.00276(6) 


-0.0050(60) 


-0.00142 


0.00069 


-0.14661 


0.00269 


43 Ca 17+ 


1.04691 


0.00298(7) 


0.00120(60) 


-0.00156 


0.00073 


-0.13908 


0.00244 


45 Sc 18+ 


1.05191 


0.00322(7) 


0.00092(46) 


-0.00171 


0.00078 


-0.13316 


0.00223 


47jjl9+ 


1.05719 


0.00347(7) 


0.00161(81) 


-0.00186 


0.00083 


-0.12769 


0.00205 


49 T ;19+ 


1.05719 


0.00345(7) 


0.00138(69) 


-0.00186 


0.00083 


-0.12772 


0.00205 


51y20+ 


1.06277 


0.00369(8) 


0.00108(54) 


-0.00201 


0.00089 


-0.12288 


0.00190 


53Q.21+ 


1.06864 


0.00398(8) 


0.00151(75) 


-0.00216 


0.00094 


-0.11840 


0.00176 


55 Mn 22+ 


1.07481 


0.00427(8) 


0.00110(55) 


-0.00232 


0.00099 


-0.11440 


0.00164 


57p e 23+ 


1.08130 


0.00457(9) 


0.00282(141) 


-0.00248 


0.00105 


-0.11050 


0.00154 


59 Co 24+ 


1.08811 


0.00487(9) 


0.00134(67) 


-0.00264 


0.00111 


-0.10728 


0.00144 


61 Ni 25+ 


1.09524 


0.00519(9) 


0.00193(96) 


-0.00281 


0.00117 


-0.10410 


0.00136 



TABLE VI: Comparison of the ground-state hyperfine splitting AE^ of Li-like ions between different 
theoretical calculations and experimental data. The values of the energies are given in meV 



Ion 




this work 


[12] 


[13J 


Exp. [57] 


14 N 4+ 


0.40376 


0.017532(8) 


0.017532(10) 


0.017667 




27^10+ 


3.6415 


1.0283(3) 


1.0281(6) 


1.0326 




45 Sc 18+ 


4.7565 


6.0631(32) 


6.063(6) 


6.0767 


6.20(8) 


57p e 23+ 


0.090623 


0.40345(64) 


0.4036(7) 


0.40470 





seen from Table IVIIIl the values of the correction Dy^~ el are smaller than the uncertainty of 
the self-energy contribution D SE . Therefore, in the case of light B-like ions one can neglect the 
contributions of the Wichmann-Kroll terms. 

In Table HXl numerical results for the individual contributions to the hyperfine splitting in light 
B-like ions are displayed. In contrast to the hydrogen- and lithiumlike sequences, where the va- 
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TABLE VII: Individual contributions to the screened self-energy correction for the ground-state hyperfine 
splitting of the boronlike ions, in units of the function Dse- 



z 


^SE 


n vr (°) 
-^se 


n vr(l+) 
-^SE 




10 


0.000 


0.630 


-0.521 


0.109 


15 


-0.001 


0.742 


-0.608 


0.132 


20 


-0.002 


0.781 


-0.648 


0.131 


25 


-0.005 


0.788 


-0.666 


0.117 



TABLE VIII: Individual contributions to the screened vacuum-polarization correction for the ground-state 
hyperfine splitting of the boronlike ions, in units of the function Dyp. 



z 


nUe- el 
-^VP 


7-)Ue— ml 
-^VP 


n WK-el 
-^VP 




10 


0.00011 


0.00066 


-0.000000084 


0.00078 


15 


0.00059 


0.0023 


-0.00000097 


0.0029 


20 


0.0018 


0.0052 


-0.0000051 


0.0069 


25 


0.0041 


0.0096 


-0.000018 


0.0137 



lence electrons are in the s states, boronlike ions have the valence electron in the pi/ 2 state; its 
electron density vanishes at the origin. Since the nuclear structure corrections 5 and e arise from 
the nuclear region, these contributions are much smaller in the B-like ions than in correspond- 
ing H- and Li-like ions. For the uncertainty of the radiative correction we prefer a conservative 
estimation as the difference of QED corrections calculated with and without screening potential. 
This uncertainty dominates for ions in the low-Z region. Evaluation of the screened radiative 
correction within the rigorous QED approach is presently underway. For high-Z ions the total 
theoretical uncertainty is mainly determined by the frequency-dependent (QED) contribution in 
the higher-order interelectronic-interaction term C(Z,aZ)/Z 2 , which is estimated to be of the 
order (aZ) 3 C(Z, aZ)/Z 2 . In Table |X] the total theoretical values of the energies AE (2pi /a) and 
wavelengths \( 2pi /^ are reported. The recoil correction is accounted for by a factor (1 + m/M)~ 3 
in Eq. © with 100% uncertainty, caused by uncalculated specific mass shift and spin-orbit recoil 
corrections [56]. 
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TABLE IX: Individual contributions to the ground-state hyperfine splitting of the boronlike ions. The values 
of the nuclear parameters are the same as in Table HI 



Ion 


A(aZ) 


8 


£ 


^rad 


B(aZ)/Z 


C(Z, aZ)/Z 2 


14 N 2+ 


1.00513 


0.000001 


0.000000(1) 


0.00017(37) 


-0.82166 


0.15006(2) 


15 N 2+ 


1.00513 


0.000001 


0.000002(2) 


0.00017(37) 


-0.82166 


0.15006(2) 


17 3+ 


1.00671 


0.000002 


0.000001 


0.00020(33) 


-0.72061 


0.11479(2) 


19p4+ 


1.00850 


0.000002 


0.000001(1) 


0.00023(29) 


-0.64222 


0.09083(3) 


21 Ne 5+ 


1.01052 


0.000003 


0.000002(1) 


0.00025(26) 


-0.57969 


0.07381(3) 


23 Na 6+ 


1.01275 


0.000005 


0.000002(1) 


0.00027(24) 


-0.52869 


0.06128(3) 


25 Mg 7+ 


1.01520 


0.000006 


0.000003(1) 


0.00029(21) 


-0.48635 


0.05179(3) 


27 A1 8+ 


1.01788 


0.000008 


0.000003(1) 


0.00030(19) 


-0.45068 


0.04443(4) 


29 si 9+ 


1.02078 


0.000010 


0.000004(2) 


0.00031(17) 


-0.42023 


0.03860(4) 


31pl0+ 


1.02392 


0.000013 


0.000005(3) 


0.00031(16) 


-0.39398 


0.03391(4) 


33gll+ 


1.02728 


0.000016 


0.000009(5) 


0.00032(14) 


-0.37113 


0.03008(5) 


35q12+ 


1.03089 


0.000021(1) 


-0.000003(11) 


0.00032(13) 


-0.35110 


0.02691(5) 


37 C1 12+ 


1.03089 


0.000021(1) 


-0.000006(14) 


0.00032(13) 


-0.35110 


0.02691(5) 


39 K 14+ 


1.03882 


0.000030(1) 


-0.000026(38) 


0.00032(10) 


-0.31768 


0.02203(6) 


41 K 14+ 


1.03882 


0.000030(1) 


-0.000062(74) 


0.00032(10) 


-0.31770 


0.02203(6) 


43 Ca 15+ 


1.04316 


0.000036(1) 


0.000016(8) 


0.00032(9) 


-0.30363 


0.02012(6) 


45 Sc 16+ 


1.04775 


0.000043(1) 


0.000014(7) 


0.00032(8) 


-0.29103 


0.01849(7) 


47 T jl7+ 


1.05260 


0.000051(1) 


0.000026(13) 


0.00032(7) 


-0.27967 


0.01708(7) 


49j;17+ 


1.05260 


0.000051(1) 


0.000023(11) 


0.00032(7) 


-0.27967 


0.01708(7) 


51yl8+ 


1.05772 


0.000060(1) 


0.000019(10) 


0.00031(6) 


-0.26941 


0.01585(7) 




1.06311 


0.000070(1) 


0.000029(15) 


0.00031(5) 


-0.26010 


0.01477(8) 


55 Mn 20+ 


1.06877 


0.000082(2) 


0.000023(12) 


0.00030(5) 


-0.25163 


0.01383(8) 


57p e 21+ 


1.07472 


0.000095(2) 


0.000064(32) 


0.00030(4) 


-0.24390(1) 


0.01300(9) 


59 Co 22+ 


1.08096 


0.000109(2) 


0.000033(16) 


0.00029(3) 


-0.23685(1) 


0.01226(9) 


61 N j23+ 


1.08749 


0.000125(2) 


0.000051(25) 


0.00029(2) 


-0.23039(1) 


0.01160(10) 



IV. DISCUSSION 

Predictions for the total energies AE'w and wavelengths A*") of the transitions between the 
ground-state hyperfine splitting components of the light H-, Li-, and B-like ions are given in 
Table |X] Due to the discrepancies in the experimental data for nuclear magnetic moments \i for 
some ions we have evaluated transition energies and wavelengths for all values of fi reported in 
Ref. [54]. The values of the nuclear spin and parity, and the root-mean-square radii are the same 
as in Table HI In the parentheses the uncertainty of the presented results is indicated. For Is and 
2s states it is mainly due to the Bohr-Weisskopf effect and must generally be considered as the 
order of magnitude of the expected error bar. For some ions we give also a second value for the 
uncertainty, which corresponds to the uncertainty of the nuclear magnetic moment. The values 



17 



TABLE X: The energies AE^ (meV) and wavelengths A^ a ^ (cm) of the transitions between the ground- 



state hyperfine splitting of the H-, Li-, and B-like ions. 



Nucleus 




AE< ls > 


A(is) 


AB< 2s > 


A< 2 ») 


A£( 2 Pi/2) 


\( 2 Pl/2) 


14 N 


0.40376 


0.21936(6) 


0.56521(16) 


0.017532(8) 


7.0720(31) 


0.0030389(34) 


40.799(46) 


15 N 


-0.28319 


0.20490(18) 


0.60510(53) 


0.016376(23) 


7.5711(105) 


0.0028419(32) 


43.627(49) 


17 


-1.8938(1) 


1.2294(2)(1) 


O.10O85(2)(l) 


0.10497(3)(1) 


1.181 1(3)(1) 


0.020459(17)(1) 


6.0602(50)(3) 


19p 


2.6289 


4.0541(7) 


0.030582(6) 


0.36364(9) 


0.34095(9) 


0.076848(49) 


1.6134(10) 


21 Ne 


-0.66180(1) 


0.93432(27)(1) 


0.13270(4) 


0.087074(34)(1) 


1.4239(6) 


0.019531(10) 


6.3481(33X1) 


23 Na 


2.2175 


4.1738(8) 


0.029705(5) 


0.40110(9) 


0.30911(7) 


0.094130(42) 


1.3172(6) 




2.2177 


4.1742(8) 


0.029703(5) 


0.40113(9) 


0.30909(7) 


0.094138(42) 


1.3170(6) 


25 Mg 


-0.85545(8) 


1.8840(5)(2) 


0.065809(19)(6) 


0.18567(7)(2) 


0.66777(25)(6) 


0.045176(17)(4) 


2.7445(10X3) 


27 A1 


3.6415 


10.215(2) 


0.012137(3) 


1.0283(3) 


0.12058(4) 


0.25755(8) 


0.48140(16) 


29 Si 


-0.55529(3) 


3.2484(10)(2) 


0.038168(12)(2) 


0.33293(13)(2) 


0.37240(15)(2) 


0.08543 1(24)(5) 


1 .4513(4)(1) 


31p 


1.1316 


8.1582(29) 


0.015198(5) 


0.84933(36) 


0.14598(6) 


0.22240(6) 


0.55748(14) 


33 s 


0.64382 


3.7624(20) 


0.032954(18) 


0.39711(25) 


0.31221(20) 


0.10583(2) 


1.1716(3) 


35 C1 


0.82187 


5.7821(64) 


0.021443(24) 


0.61780(81) 


0.20069(26) 


0.16687(3) 


0.74298(15) 


37 C1 


0.68412 


4.8143(67) 


0.025753(36) 


0.51439(85) 


0.24103(40) 


0.13891(3) 


0.89258(18) 


39 K 


0.39147 


3.873(12) 


0.03202(10) 


0.4225(15) 


0.2934(11) 


0.11670(2) 


1.0625(2) 




0.39151 


3.873(12) 


0.03201(10) 


0.4226(15) 


0.2934(11) 


0.11671(2) 


1.0624(2) 


41 K 


0.21487 


2.132(13) 


0.05816(35) 


0.2326(16) 


0.5331(37) 


0.064053(12) 


1.9356(4) 




0.21489 


2.132(13) 


0.05816(35) 


0.2326(16) 


0.5330(37) 


0.064059(12) 


1.9355(4) 


43 Ca 


-1.3176 


13.025(8) 


0.0095191(57) 


1.4340(10) 


0.086460(60) 


0.40134(6) 


0.30892(5) 


45 Sc 


4.7565 


54.613(25) 


0.0022702(11) 


6.0631(32) 


0.020449(11) 


1.7115(2) 


0.072440(10) 


47 Ti 


-0.78848(1) 


10.957(9) 


0.011316(9) 


1.2259(11) 


0.10114(9) 


0.34905(5) 


0.35520(5) 


49 Ti 


-1.1042 


14.617(10) 


0.0084821(59) 


1.6355(13) 


0.075810(60) 


0.46554(6) 


0.26632(3) 


51 V 


5.1487 


78.168(42) 


0.0015861(9) 


8.8099(54) 


0.014073(9) 


2.5248(3) 


0.049107(6) 


53 Cr 


-0.47454(3) 


9.5797(72)(6) 


0.012942(10)(1) 


1.0871(9)(1) 


0.11406(10)(1) 


0.31367(4)(2) 


0.39527(5)(2) 


55 Mn 


3.4687 


71.525(39) 


0.0017334(9) 


8.1688(50) 


0.015178(9) 


2.3699(3) 


0.052315(6) 




3.4532(13) 


71.206(39X27) 


0.0017412(9)(7) 


8.1323(50)(31) 


0.015246(9)(6) 


2.3594(3)(9) 


0.052550(6)(20) 


57 Fe 


0.090623 


3.5109(49) 


0.035314(50) 


0.40345(64) 


0.30731(49) 


0.11787(1) 


1.0519(1) 




0.090764 


3.5164(49) 


0.035259(50) 


0.40408(64) 


0.30683(48) 


0.11805(1) 


1.0503(1) 




0.09044(7) 


3.5038(49)(27) 


0.035386(50)(27) 


0.40264(64)(31) 


0.30793(49)(24) 


0.11763(1)(9) 


1.0540(1)(8) 


59 Co 


4.627(9) 


115.35(8)(22) 


0.0010749(7)(21) 


13.334(10)(26) 


0.0092987(69)(181) 


3.9084(5)(76) 


0.031722(4)(62) 


61 Ni 


-0.75002(4) 


24.418(23X1) 


0.0050777(49)(3) 


2.8385(30)(2) 


0.043680(47)(2) 


0.83617(10)(4) 


0.14828(2)(1) 



for 



ithiumlike 45 Sc 18+ and boronlike 45 Sc 16+ , 57 Fe 21+ ions coincide with our previous results 



44, 



45|]. Due to the lack of experimental data one can not make a detailed comparison for the 



ions under consideration. 

Table IXl shows excellent accuracy for the values of the hyperfine splitting of H-, Li- and B-like 
ions. For the H- and Li-like ions the limitation of the total accuracy is set by the Bohr-Weisskopf 
correction. However, this uncertainty can be considerably reduced in the specific difference of the 
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0] 



ground-state hyperfine structure values of H- and Li-like ions with the same nucleus 

A'E = AE {2s) - £A£ (ls) . (25) 

The parameter £ has to be chosen to cancel the Bohr-Weisskopf correction 

1 A™{*Z){l-8™) 
t = 8AM( a Z)(l-6M) f(aZ) - (26) 

The function f(aZ) is defined by the ratio of the Bohr-Weisskopf corrections 

£ (2s) 

fiflZ) = . (27) 



This function can be calculated to a rather high accuracy, because it is determined mainly by the 
behavior of the wavefunctions at the atomic scale and thus almost independent of the nuclear 
structure [9]. In Table IXJ we present the numerical results for the parameter £ and for the specific 
difference A'E between the hyperfine splitting values of H- and Li-like ions. The theoretical 
accuracy of the presented values A'E is better than 0.01% and the given uncertainty is determined 
by the limited knowledge of the nuclear magnetic moments. 

Let us summarize: Ab initio QED calculations of the ground-state hyperfine splitting of H-, Li-, 
and B-like ions in the middle-Z region have been performed. The evaluation incorporates results 
based on a rigorous treatment of first-order many-electron QED effects and on the large-scale 
CI-DFS calculations of the second- and higher-order electron-correlation effects. The one-loop 
radiative corrections have been evaluated to all orders in aZ. The screening QED effect in Li- 
and B-like sequences have been taken into account utilizing a local Kohn-Sham potential. The 
Bohr-Weisskopf correction has been calculated employing the single-particle nuclear model. As 
the result, the most accurate values for the hyperfine splitting of H-, Li-, and B-like ions under 
consideration have been obtained. The specific difference of the hyperfine splitting values of 
H- and Li-like ions, where the uncertainties associated with the nuclear- structure corrections are 
significantly canceled, has been evaluated with an accuracy better than 0.01%. 
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TABLE XI: The parameter £ and specific difference A' ' E (meV) between the hyperfine structure values of 
H- and Li-like ions, defined by Eqs. (l26l ) and (|25T ). respectively 



Nucleus n/fJ-N £, A' E 



14 N 


0.40376 


0. 


12531 


-0.009955 


15 N 


-0.28319 


0. 


12530 


-0.009298 


17 o 


-1.8938(1) 


0. 


12539 


-0.04919 


19p 


2.6289 


0. 


12550 


-0.1451 


21 Ne 


-0.66180(1) 


0. 


12561 


-0.03029 


23 Na 


2.2175 


0. 


12574 


-0.1237 




2.2177 


0. 


12574 


-0.1237 


25 Mg 


-0.85545(8) 


0. 


12588 


-0.05148 


27 Al 


3.6415 


0. 


12603 


-0.2592 


29 Si 


-0.55529(3) 


0. 


12619 


-0.07699 


31p 


1.1316 


0. 


12637 


-0.1816 


33 S 


0.64382 


0. 


12656 


-0.07904 


35 C1 


0.82187 


0. 


12676 


-0.1151 


37 C1 


0.68412 


0. 


12676 


-0.09587 


39 K 


0.39147 


0. 


12720 


-0.07006 




0.39151 


0. 


12720 


-0.07007 


41 K 


0.21487 


0. 


12720 


-0.03856 




0.21489 


0. 


12720 


-0.03857 


43 Ca 


-1.3176 


0. 


12743 


-0.2258 


45 Sc 


4.7565 


0. 


12768 


-0.9098 


47 Ti 


-0.78848(1) 


0. 


12794 


-0.1759 


49 Ti 


-1.1042 


0. 


12794 


-0.2347 


51y 


5.1487 


0. 


12822 


-1.213 


53 Cr 


-0.47454(3) 


0. 


12851 


-0.1440 


55 Mn 


3.4687 


0. 


12881 


-1.044 




3.4532(13) 


0. 


12881 


-1.039 


57 Fe 


0.090623 


0. 


12912 


-0.04989 




0.090764 


0. 


12912 


-0.04996 




0.09044(7) 


0. 


12912 


-0.04979(4) 


59 Co 


4.627(9) 


0. 


12945 


-1.599(3) 


61 Ni 


-0.75002(4) 


0. 


12979 


-0.3307 
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